EXPLICIT INVERSIONS OF CERTAIN MATRICES II 



RUIMING ZHANG 



Abstract. In this note, we apply kernel polynomials to find the explicit in- 
verses for some some Hankel matrices associated with g-orthogonal polynomi- 
als. 



1. Introduction 

In the theory of orthogonal polynomials, We could calculate the determinants 
of some Hankel matrices once we know the three term recurrence relation for the 
associated orthogonal polynomials and vice versa. It is well-known that the kernel 
polynomials of the orthogonal polynomials encode important information about 
the associated Hankel matrices. These matrices are generalizations of the Hilbert 
matrices. In this note we invert some Hankel matrices associated with some q- 
polynomials by using the kernel polynomials. 

The rest of the introduction could be found in [4] , we copy it here for the conve- 
nience of the reader. 

Theorem 1.1. Given a probability measure [i onR with a support of infinite many 
points. Let us consider the Hilbert space of [i-measurable functions 



X:=\f{x)\ / \f(x)\ 2 dn(x) < 



(1.1) 

with the inner product defined as 

(1-2) (f,g) := I f(x)gjx)dfi(x), f,geX. 



Assume that {w n (x)}'^' =Q is a sequence of linearly independent functions in X with 
wq(x) = 1. Let 



(1.3) 



(1.4) 



Wj(x)wk(x)dfi(x), j,k = 0, 1. 



n„ := 



aio an 
\ a„o a n \ 



OL\ n 



n e NU{0}. 
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(1.5) 



A„:=detn„, n G N U {0} . 



Then the n-th orthonormal function with positive coefficient in w n {x) is given by 
the formula 



( «00 "01 



(1.6) Pn(x) 



: (let 



&02 



a w an «12 



a„ a„i a n2 
\ w (x) wi(x) w 2 {x) 

for n£tf, with 

(1.7) p (x) = w (x) = 1. 

Furthermore, the coefficient of p n (x) in w n (x) is 



Otln 



w n (x) J 



(1.8) 



7n : = 



A„_i 



Proof. The proof is the same as for the case w n (x) = x n , which could be found in 
any orthogonal polynomials textbooks such as [3]. □ 

Corollary 1.2. For n £ N ; we have 

n 



(1.9) 



k=l 



In 



Proof. This is a trivial consequence of (|1.7p and (|1.8j) . 
Lemma 1.3. Let 



□ 



(1.10) 



k n (x,y) := ^2pk{x)pk(y), neNU{0}. 



Then, for any it(x) in the linear span of {wk{x)}^ , we have 



(1.11) J ir(x)k n (x,y)dfj.(x) = n(y). 

Proof. To see (|l.lip . just expand n(x) in Pk(x), k = 0, 1 . . . , n. □ 
Lemma 1.4. For each n S NU{0} ; the function k n {x, y) satisfying is unique. 

Proof. Suppose there are two such functions h n (x,y) and k n (x,y) , then, 

(1.12) 0< |M-,y)-M-,y)|| 2 

=(h n (; y) - kn(;y),h n (;y) - k„(-,y)) 

=(h n {-,y) - k n (-,y),h n (-,y)) - (h n (-,y) - k n (- , y) , k n (- , y)) 
=0, 

which is a contradiction. □ 
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Lemma 1.5. Let 

(1.13) (fok)o<j,k<n ■■ = n" 1 , neNU{0}. 
Then, 

n 

(1.14) k n (x,y) = Pj k Wj(y)w k (x). 
Proof. Let 

n 

(1.15) f(x) = }^UkWk(x), 

k=0 

then, 

n 

(i-i6) (/(•), E PjkwffiM-)) 

j,k=0 

n 

= E W ™K(').K-.!/)) 

m=0 
n n 

= E Um E l3jk w j(y)( w m,W k ) 
rn=0 j,fc=0 
n n n 

= E Um E E Pjk a km 

m=0 j=0 k=0 



By Lemma [Ql we have 

(1.17) k n (x,y) = E PjkWj(y)w k (x). 

j,k=0 

Corollary 1.6. The kernel in Lemma f"OI is also given by 



(1.18) k n {x,y) = -— det 



/ 1 u>i(y) • 

1 aoo «oi 

I0i(:r) a w an 

\ w n (x) a n0 a n i 



forn 6 NU {0}. 
Proof. Since 

n 

(1.19) k n (x,y) = E P]kWj{y)wk{x) 

j,k=0 

with 
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Then, 
(1.21) 



n„(fc,j) n„(fc,j) 



Jk detn„ A„ 
where n n (fc,j) is the (k,j)-th co-factor. Therefore, 



(1.22) 

It is clear that 
(1.23) 



I « 

k n {x,y) = — ^2 u n(k,j)wj{y)w k (x). 



j,k=0 



^2 Tl n(k,j)w j (y)w k (x) 
j,k=o 



(W(y)) 

w(x) n„ 



by direct determination expansion, which is 



(1.24) 
where 

(1.25) 

and 
(1.26) 



k n (x,y) 



1 



W(x) 



(W(y)) 

w(x) n„ 



( 1 \ 

\ w n (x) J 



(W(y)) = ( 1, Wl (y) w n (y) ) 



□ 



Lemma \l .51 enables us to compute the inverse the Gram matrix in terms of the 
orthonormal functions {pn(x)} n — - 

Corollary 1.7. Assume that {w n (x)}^ =Q , {p n {x)}^ =0 andll n = {a 3 k)o<jM<n as 
in Theorem \l.l[ Suppose we have two families of linear junctionals {wfc}^L an d 
{ v k}kLo over li near space generated by {w n (x)}'^' =0 with 



(1.27) 
and 



Uj(w k ) = Sjk, 



(1.28) 

for j, k = 0, 1, . 



Vj(w k ) = S jk 



Then, 



(1.29) 

where 
(1.30) 



/%* = U k(Pm{x))Vj(p m (y)), 



ifljk) 



'jk)0<j,k<n 



n: 
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Proof. From Lemma 11.51 we have 

n n 

(1.31) ^2 PjkWj(y)w k {x) = 2J Pm{y)p m (x). 

j,k=0 m=0 

Then we apply the functional Uj and Vk both sides of the above equation, the claim 
of the corollary follows. □ 



2. Main Results 



As usual, [HE], 



(2.1) 

and 
(2.2) 



(a;g)oc := Y[(l-aq m ), a e C, q 6(0,1), 



m=0 



(a;g) r 



(a;g)c 



(ag m ;g)oo ' 

We will use the following short hand notations 



m 6 Z. 



(2.3) 
and 



m 
j 



(<?;<?)r 



j <m, j, m 6 N U {0} , 



(2.4) (ai,a 2 , ...,a„;q) m : = _Q (a^; g) TO > n 6 N, to 6 Z, ai,...,a n 6 



fc=i 



The basic hypergeometric function r s with complex parameters a\, a r ; &i, b s 
is formally defined as 



(2.5) 



"i . ...\. (ai,—,a r ;q) m ( { („-i), l /2 N i s+1 r 



E 



„^ fo & i,-> &«;?)» 



(-i)V 



6i, ...,6 S 

The g-binomial theorem is the following summation formula 
(2.6) 



and the Ramanujan summation formula is 
(2.7) 



n— — oo 



(6;g)„ (b,b/az,q/a,z;q) 
Theorem 2.1. For neNU {0} ; the matrix 

(2.8) ^ (a<?;9)j+fc 

has determinant 



< z < 1. 



{abq 2 ;q)j +k J fe=0 
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(2.9) 



det 



(abq 2 ;q) ]+k J Jk=0 



-a 



n(n+l)/2 n(n+l)(2ra+l)/6 



n 

k=l 



{q,aq, bq,abq;q) k 
(abq, abq 2 , abq 2 , abq 3 ; g 2 )* 



The matrix (|2 .8[) is invertible for 

(2.10) q k ,aq k ,bq k ,abq k ^ 1, fc G N, a,g^0, 

its inverse matrix (7jfc)™j. =0 nas elements 

(_l)i+feg[i(i+i)+fe(fe+l)]/2 



(2.11) 



E 



(aq;q)j(aq;q)k 

(1 - abq 2m+1 ){aq,abq;q) r 



m-max(j,fc) 



m 




m 




9 


k 



(l-abq)(q,bq- q) m (aq^ +k+1 ) m 



(abq^q^abq^q)* 



Remark 2.2. Let a = 6 = 1 in H2.8J1 . then take g — > 1 we would get the classical 
Hilbert matrix. 

Theorem 2.3. For n G N U {0} and g ^ 0, £/ie matrix 

(IT" 



(2.12) 
/ias determinant 



j+kq 



j,k=0 



(2.13) det((< Z «+ 1 ; g ) J+fc9 -^ 

TTie matrix (|2.12p is invertible for 

(2.14) -a£N, g^l.teN 
i/«e inverse matrix ('Jjk)™ k=o has elements 



I\Li(q,q a+1 ;qh 



k=0 on(n+l)(4n-l)/6 



(2.15) Tifc 



9 



[j(j-i)+fc(fe-i)] 



(« a+1 ;g)j(? a+1 ;g)fc 



E 



m— max(j,fc) 



(<z a+1 ;<?)™ 






m 




i 


9 


k 



In our last result, we use the following notations 

(2.16) x = [x\ + {x} , L^JeZ, {x}e[0, l). 

Theorem 2.4. i<br n G N U {0} and q ^ 0, the matrix 



(2.17) ( 

has the determinant 

(2.18) det 



[1 + (-!)'■+*] _ 



u+ky 



j,k=0 



q )\i+k 



IILo(g;g)fc 

7 n(n+l)(2n+l)/6 ' 



j,fc=0 
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For 

(2.19) q k ^ 1, k £ N, 

the matrix (|2.17p «s invertible and has inverse 
(2.20) 



E 



m 




m 


j 


1 


k 



(g, q 2 ) (g, q 2 ) cos 



(m-j)j (m-fc)Tr \ 



(-l)-(g;g) m g- m -(S-(5) 



V 



/ j,fc=0 



3. Proofs 

Given a polynomial f(x), we define the g-difference operator as [lj [2] 

f(x) - f(qx) 



(3.1) (D q f)(x) : = 

Let Wk{x) := x fe for fc = 0, 1, .... then we have 
(3.2) 
and 



1 - g" 

D q w n {x) = -j— — w n -i(x), 



(3.3) 



«*(/(*)) : = «*(/(*)) := lW)(*)L =0 > 



For this polynomial sequence we take 
(3.4) 

where f(x) is a polynomial in variable x. 

3.1. Proof for Theorem 12. 1L The little g-Jacobi polynomials {p n (x; a, 6|g)}^L 
are defined as [TJ [2] 

/ q-n abq n+1 \ 

(3.5) Pn{x;a, b\q) : = 2 4>1 I ' ;q;qxj, < g, ag, 6g < 1 

for n > 0, and we assume that 

(3.6) p_i(x;a, b\q) : = 0, po(z; a, b\q) := 1. 

The little g-Jacobi polynomials {p n (x; a, b\q)}™_ have the orthogonal relation 



(3.7) 



E 7773 Pm(q k ;a, b\q)p n (q k ;a, b\q) = h n (a, b\q)5„ 



k=0 

for m.n > with 



(3.8) 



fr„(a,%) : = 



{abq 2 :q)oo (1 - abq){aq) n (q 7 bq;q) n 



(ag; 17)00 (1 - a6g 2n+1 ) (ag, a&g; g)„ ' 
The little g-Jacobi polynomials satisfy the following difference relation 

/ anl n fe / tn q k (q~ n ;q) k {abq n+1 ;q) k k k 

(3.9 D K p n {x; a, b\q) = -j- 7 p n -k(x; aq K ,bq K \q), 

q (1- qr(aq;q)k 
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for n,k = 0, 1... and n > k. The orthonormal polynomial 

(3.10) p n (x) = ^)>n(x;ag,bq\q) 

^h n (a,b\q) 

have the leading coefficient 



(3.11) 7 n 



(aq; q) 00 (abq, abq 2 , abq 2 , abq 3 ; q 2 ) n 
(abq 2 ;q) QO (q, aq, bq, abq; q) n a n q n2 



for n = 0, 1.... 

The moments are given by the formula 



(3.12) Mn = 



(bq;q) m (aq) m q* 



or 

(abq n+2 ;q) c 



(3 - 13) M " = {aq^q) c 

by using the g-binomial theorem. 
For any n = 0, 1, the matrix 



(3.14) H n 



f (abq^+ 2 ;q) 



V {aqi+^-q)^ J ] k=0 

has the determinant 

(3 15) det fW +k+2 ;q)oo\ : _ (wb.r.q).^" V"-" +h ' 



\{aq3+ k + 1 ;q) 00 J ]k=0 \ (aq;q)oo 

(q, aq, bq, abq; q) k 



n+l 



X 



n 



(abq, abq 2 , abq 2 , abq 3 ; q 2 )k ' 



(3.16) det ( ^ q i q \+ k ) = a "(n+l)/2 9 n(n+l)(2n+l)/6 

\{abq 2 ;q) 0+k J ]k=0 

x tt (q,aq,bq,abq;q) k 
1J- (abq, abq 2 , abq 2 , abq 3 ; q 2 ) k 

by applying Theorem ll.il 

Under the condition < q,aq,bq < 1, the elements of H^ 1 — ((3jk)o<j,k<n are 
given by 

, al „s o q j+k (aq;q)oo ^ (1 - abq 2m+1 )(aq,abq;q) m 



(abq 2 ;q) x m= ^ (jk) (1 — abq)(aq) m (q, bq; q) r 

x (g- m ;g) J -(abg m+1 ;g) J - (g- m ;g) fc (afeg m+1 ;g) fc 



(q\q)j(aq;q)j (q;q)j(aq;q) k 
for j, fe = 0, 1, n. Then, the matrix 

(3.18) f (a9; ^ +fe 



(abq 2 ;q) J+k J ]k=Q 
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has the inverse matrix (jjk)™ k=0 with 

(_l)i+fc g b'0'+i)+fc(k+l)]/2 



(3.19) 



Ijk 



(aq; q) 3 (aq; q) k 

y> (1 - abq 2m+1 )(aq,abq;q) m 

^ (1 — abq)(q, bq; q) m (aqi +k+1 ) m 

m—max[j,k) 



m 
j 



m 
k 



(abq m+1 ;q) J (abq m+1 ;q) k . 



Since the formulas (|3. 16fl . (|3.18p and ([3.19P contain only rational functions of a, b 
and q k , k e N, the original restrictions in (|3.5p could be dropped and Theorem 12.11 
follows. 

3.2. Proof for Theorem 12. 31 The g-Laguerre polynomials { L^\x; q) \ are 

I J n=0 

defined as [IJ [2] 



(3.20) 

L^Hx-q) 



{q a+1 ;q)n {g- n ;q) k q^ 1 ){-x) k g^+^ k 



q e (0, 1), a > -1 



for n > 0, and we assume that 

(3.21) L^(x;q): = 0, L^\x; q) := 1. 

The g-Laguerre polynomials have an orthogonal relation 



(3.22) 



oo ka-\-k 



k— — c 



(-cq k ;q) c 



for m, n = 0, 1, ... with 



(3.23) 



/i„(c, a\q) : 



(g Q+1 ,-c, -c ^jgOoc (g;g)«g" ' 

The g-Laguerre polynomials satisfy the following difference relation 



(3.24) D k L^ (x; q) = ^^L^ (q k x; q) 

for n, k = 0, 1.... The orthonormal g-Laguerre polynomial 



(3.25) 

has the leading coefficient 



*n(a;) 



L ( n\x;q) 
y/h n (c,a\q) 



(3.26) 

for n = 0, 1. 



,(a+n+l/2)n 



In 



(g"+ 1 ,-c,-c- 1 g;g) c 



V(g,g Q+1 ;9)n V (9,-«z a+1 ,-c-V Q ;<7)oo 
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The moments are given by 

00 r n n k(u+n+l) 

(3.27) M n= ! 



k= 



c" 



(-cq k ;q)oo 

CO 

{-c;q) k q k{a+n ^ 



1) 



which could be summed by the formula (|2.7[1 . 

c n (q,-cq a+n+1 1 -q- a - n /c;q) c 



(3.28) fl n 



(-c;g)oo (-g/c,g a +" +1 ;g)oo 
= c"( g «+ 1 ; g) n (-g- a -"/c; g) n (g, -cg^ 1 , - g -°/c; g) 
(-cg«+ 1 ;g)„ (-c,-g/c,g"+ 1 ;g) O0 
(g" +1 ;g)» (g,-cg a+1 ,-g- a /c;g) 
qna +n(n+i)/2 (-c, -q/c, q a+1 ;q) O0 

Then for any n = 0, 1, the matrix 

(; , 2<)i //( , (<7 a+1 ;<z);+fc (g,- C g«+ 1 ,-g-«/c;g)\" 



^[(j+fc)«+0'+fe)0+fe+i)/2] (-c,-?/c,g a + 1 ;g) 00 y i)fe=0 
has the determinant 

(3.30, det e. - fyf^-'-f*) " +1 n 

or 

by applying Theorem ll.il By applying Corollary ??, we obtain the inverse matrix 
ff^ 1 = (&*)£*=<> with 

(g Q+1 , -c, -cT 1 ?; g)oo g Q (J'+ fe )+J 2+fc2 



(3.32) ifc 



(g, -cg Q+1 ,-c x g ";g)oo (g;g)j(g;g)fe 

y (gj g) m g m (g" +J+1 ; qhn-j (q a+k+1 ; q) m -k 



then, the inverse matrix (jjk)™ k=o °f 



(3.33) ((g a + 1 ;g)^,fl-( i t fe ) > l" 



has element 

_ gi'-J+fc'-fc " (g;g) m g m (g a + 1 ;g) m - J (g a+1 ;g) m - fc 

1 J 7jfc (g;g),(g;g) fc ^, a (f +1 ;«) m (g;g) m _,(g;g) m _ fc 

g b'0'-i)+fc(fc-i)] " 



(g^gMg^g)* „ ( . fe) 



te a+1 ;g)mg m 


m 




m 


(g;g)m 


j 




k 
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or 



(3.35) t,-, 



q 



b'0'-i)+Jfc(fe-i)] 



(q a+1 ;q)j(q a+1 ;q)k 



E 



m— max(j.fc) 



(q a+1 ;q) m q m 


m 




m 


(q;q) m 


i 


Q 


k 



Since the formulas p,3ip . Q3.33P and Q3.35P only involve rational functions of q a 
and q k , k £ N, thus the restriction q may be dropped to get Theorem 12,31 



3.3. Proof for Theorem 12. 4L The discrete g-Hermite polynomials II are defined 
as [2] 

/ n —n „— n+1 „2 

(3.36) h n (x;q) : = x n 2 cb n ( q 'J ^ 
for n > and we assume that 

(3.37) ft_i(a:;g): = 0, h {x; q) := 1. 
The polynomials {/j n (a;;g)| satisfy the difference equation 

I J n=0 

(3.38) 



■ - (q- q) n q- jn+j{j+ls >/ 2 ~ 
iy a -iK{x\q) = —j-, Zvf7I7Z\ K-j{x\q), 



(i - q) j (q;q)n- 3 

They also satisfy the orthogonality 
(3.39) 

oo 

\h m (cq k ]q)h n (cq k \q) + h m (~cq k \q)h n (-cq k -.q) w(cq k )q k = h n (c\q)S n 

k— — oo 

for m,n > 0, where 

(3.40) w{x) 1 
and 

(3.41) h n {c\q) : 



{-x 2 ;q 2 )oo ' 
(q 2 ,-c 2 q,-c- 2 q;q 2 ) oc (q;q) T 



(q, -c 2 ,-c 2 q 2 ;q 2 )oo q n2 
for some c > 0. thus the orthonormal polynomial 

h„(x;q) 



(3.42) 

has the leading coefficient 
(3.43) 



h n (x) : 



y/h n {c\q) 



y/h n {c\q) 
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The moments of the measure could be calculated via formula [22 



(3.44) 



Mn = Yl [(cq k ) n + (-cq k T]w(cq k )q k 
k— — oc 

q (n+l)k 



[1 + (-!)»] c » 



k— — c 



(—c 2 q 2k ;q 2 ) 



[l + (-l)"]c" 
(-c 2 ;g 2 )oo 



E (~ c W) 



A: '7 



(n+l)fe 



A— — oo 

_ [1 + (-1)"] e" (g 2 , -c- 2 g-" +1 , -c 2 q n+1 ;q 2 ) 

(-cW)oo (-^V, g 2 )oo 

_ [1 + (-1)"] g -" 2 / 4 ( g 2 , -c- 2 g, -c 2 g; g 2 )^ 
(-c 2 , -~c~ 2 q 2 , g n+1 ; q ,2 )oo 
Thus, for n = 0, 1, the matrix 

' 1 + (-lp' +fe ] q- ( ^ 2 ^(q 2 , -c- 2 q, -c 2 q; q 2 ) c 



(3.45) H n 
has determinant 



(- c 2 ,- c -V,^' +fc+1 ;9 2 )c 



j,fe=0 



(3.46) 



dot H„ 



IILo(g;g)fc 

7 n(n+l)(2ri+l)/6 



, (g 2 ,-c 2 g, -c 2 g;g 2 ) 
' (q,-c 2 ,-c- 2 q 2 ;q 2 ) c 



n+l 



Therefore, 
(3.47) det 



[l + (-l) J+fe ] - 



IILo(g;g)fc 

7 n(n+l)(2n+l)/6 ' 



j,fc=0 

For n = 0, 1, the matrix H^ 1 = (/3j&)™^ =0 nas element 

(3.48) 0jk = ( -iy+y 2 +* 2 JvV c ~ 2g2 / g2 2 r 

2(g 2 , -c*q,-c z q;q z )oo 
q^-U+^hm-Mqyhm-k&q) 



n 

x E 

m-max(j,fc) 



(<z;g)r 



m 
3 



m 
k 



From the generating function 

(-xt; q) c 



(3.49) 
we get 
(3.50) 

Then, the matrix 
(3.51) 



q^)t n ~ 



(-t 2 ;q 2 ) c 



h n {Q\q) 



E 7 — \-hn{x;q), 



(g;g 2 )oocos^ 
(<Z 1+ ";g 2 )oo^' 



i,fe=0 
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has inverse matrix 
(3.52) 

/ 

n 

E 

m-max(j.fc) 

V 



m 




m 


3 




k 



(q, q 2 ) (q, q 2 ) cos (m 2 j)7r cos 



(-l) ro (g;«) m g- m -®-(") 



(m — fc)7r \ 
2 

/ 



j,k=0 



Since formulas <|3.47|) . I|3.5ip and p.52p contain only rational functions of q k , k E N, 
thus the restriction may be dropped and Theorem 12.41 follows. 
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